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The spin echo approach is extended to include bio-complexes for which the interaction with dy¬ 
namical noise is strong. Significant restoration of the free induction decay signal due to homogeneous 
(decoherence) and inhomogeneous (dephasing) broadening is demonstrated analytically and numer¬ 
ically, for both an individual dimer of interacting chlorophylls and for an ensemble of dimers. This 
approach is based on an exact and closed system of ordinary differential equations that can be easily 
solved for a wide range of parameters that are relevant for bio-applications. 
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The femtosecond spectroscopic and measur¬ 
ing technologies, developed during last few 
decades, have allowed experimental scientists 
to examine very rapid dynamical processes, in¬ 
cluding those in biological systems. In partic¬ 
ular, it was demonstrated that, even at room 
temperature, photosynthetic bio-complexes ex¬ 
hibit collective quantum coherence (CQC) dur¬ 
ing primary electron transfer (ET) processes 
that occur on the time-scale of some hundreds 
of femtoseconds T]. 

The CQC is resulted from the fact that 
the primary processes of exciton transfer and 
charge separation are so rapid (on a time-scale 
of a few picoseconds) that the environment does 
not have time to recombine the exciton and de¬ 
stroy the CQC. 

One of the methods for studying the CQC 
effects in bio-systems is the spin echo spec¬ 
troscopy, which was initially developed in the 
nuclear magnetic resonance This spin echo 
procedure allows one to reduce the effects of 
the “inhomogeneous broadening” (dephasing) 
and to increase the time of the signal gener¬ 
ated by the transverse magnetization (free in¬ 
duction, FI). Recently, the spin echo technique 
has also been used to analyze quantum co¬ 
herent effects in an ensemble of effective 1/2 
spins (chlorophyll-based dimers) in photosyn¬ 
thetic bio-complexes. (See [§], and references 
therein.) 

During the last decade, it was recognized 
that the spin echo technique can be success¬ 
fully applied not only to an ensemble of spins, 
but also to an individual spin, or to an individ¬ 


ual two-level quantum system (TLS) Q. In this 
case, inhomogeneous broadening is absent. But 
still one can partly suppress the broadening of 
the FI decay resulted from a pure dynamical 
(time-dependent) noise. In [4|-|7| the spin echo 
approach was applied to a single quantum two- 
level system - the superconducting qubit. In 
this case, the qubit is considered in the so-called 
diagonal representation, with the main char¬ 
acteristic being the energy gap (usually tens 
of gigahertz) between the ground and the ex¬ 
cited states. The dynamical noise, acting on 
the qubit, is generated by the time-fluctuations 
of the electron charge, bias current, external 
magnetic flux, and other sources, and these all 
are relatively weak. 

In contrast to a superconducting qubit, an in¬ 
dividual effective two-level quantum system (a 
dimer) in a bio-complex is usually characterized 
in the so-called site representation, using such 
main parameters as the energy gap (redox po¬ 
tential) between the excited states of the chloro¬ 
phylls realizing this dimer, and the matrix ele¬ 
ment of the dipole-dipole or the exchange inter¬ 
actions between these two excited states. Also, 
the dimer usually experiences a strong interac¬ 
tion with the protein environment, caused by 
dynamical noise (characterized by the recon¬ 
struction energy), that must be taken into ac¬ 
count in bio-applications of the spin echo tech¬ 
nique. 

In this Letter, we analyze analytically and 
numerically the spin echo technique for both an 
individual dimer and for an ensemble of dimers 
in bio-complexes, for the case of strong inter- 


2 


action with the dynamical noise. We show 
how both the dynamical and the inhomoge¬ 
neous broadening can be suppressed by the spin 
echo pulses. Our conclusion is that even for 
strong dynamical noise, the spin echo approach 
can serve as the useful complementary spec¬ 
troscopic technique for characterizing the bio¬ 
complexes that include both individual dimers 
and an ensemble of dimers. 

We consider a dimer composed of the excited 
states of two chlorophylls, Chll and Chl 2 . We 
assume that each chlorophyll experiences a di¬ 
agonal noise, provided by the protein environ¬ 
ment, which is described by the random vari¬ 
able, £(i). In the site representation, the Hamil¬ 
tonian of the system can be written as follows: 
H = £i|l)(l| + £ 2 | 2 )( 2 | + (1/2) (V 12 11) <2| + 
/i.c.)+£(f)(Ai|l)(l| + A2|2)(2|). We also assume 
that noise is produced by the stationary ran¬ 
dom telegraph process (RTP) with: (£(f)) = 0 , 
(£W£(i')>, = X(t - t'), where, y(f - t') = 
cr 2 e _2 ' T I* — * I, is the correlation function; ct, 2y, 
and Aj., 2, are the amplitude of noise, the decay 
rate of the correlation function, and the inter¬ 
action constants with noise, correspondingly. 

In the diagonal representation of the unper¬ 
turbed Hamiltonian, we obtain the total Hamil¬ 
tonian for the effective TLS, 

T~L = —I + -Uct, + —D\ z t;(t)<j z + 

^£>A,±£(i)(cos ipa x + sin ipo y ), (1) 

where, <J x , y ,z, are the Pauli matrices, Ao = 
Ei + E2 + (Ai + A2)£(f), A = Ai — A2, H = 
\/(Ei - E 2 ) 2 + \V 12 \ 2 , D x , z = A cos 9 , and 
D \ y _l = Asind. We set: cos 9 = {Ei — E 2 )/fl. 

The dynamics of a TLS is described by two 
rates: the longitudinal relaxation rate, Ti = 
Ti , and the transverse relaxation rate, T2 = 
T 2 _1 . When the noise is weak, and the con¬ 
dition, t c <C Ti,T 2 , is satisfied (where r c = 
1/(27) is the correlation time of the noise fluc¬ 
tuations), one can apply Bloch-Redfield (BR) 
theory {§:, @|. In BR theory, the transverse re¬ 
laxation rate, T2 = T^, + Ti/2, where T v is 
the so-called “dephasing” rate. In terms of 
the spectral density of noise, S(w), these rates 
are defined as follows /$: Ti = nD 2 ± S(fl), 
r^, = nD\ z S( 0 ). Using the spectral density 
of RTP, S(ui) = 27<7 2 /7r(47 2 + w 2 ), we obtain 
the relaxation and dephasing rates provided 
by BR theory: Ti = 27U 2 sin 2 9 /{ 4 rf 2 + SI 2 ), 
r tp = (y 2 / 2 "/) cos 2 9 , where the renormalized 
interaction constant, v = Act, is introduced. 


To study the quantum decoherence, we 
present the density matrix as, p(t) = (I + n (t) ■ 
er)/ 2 , where, n (t) = Tr(p(t)er), is the Bloch 
vector. Instead of the Liouville-von Neumann 
equation for the density matrix, ihp = \H , p], it 
is convenient to employ the equation of motion 
for the Bloch vector (we set h= 1 ): 

j=llx n + (£(t)/cr)u) x n. (2) 


Here ui = v(sin 9 cosip, sin 0 simp, cos 6 ), and 

n = (0,0,H). 

Using the differential formula for the RTP 


(I + 2 7 )(£(W;£(t)]) = 

( 3 ) 


where, R[t;^(r)], is an arbitrary functional, we 
obtain from Eq. © the closed system of differ¬ 
ential equations: 

- =0 x (n) + uj x (n«), ( 4 ) 

=& x ( n? > + w x (n) - 27(1/), ( 5 ) 

at 

where, (n ? ) = {£(t)n)/a. The average, (...), is 
taken over the RTP. 

Homogeneous broadening. To characterize a 
dimer, we introduce the dimensionless parame¬ 
ter, e = | tan 0 | = \Vi 2 /(Ei—E 2 )\. When e<l, 
we will call the dimer “weakly coupled”. In 
the opposite case, e > 1, the dimer is called 
“strongly coupled”. 

We first consider a weakly coupled dimer, so 
one can neglect the effects of relaxation. Intro¬ 
ducing the complex vectors: (ro(t)) = (n x (t)} + 
i(n y (t)) and (m^(t)) = (ra£(t)) + i(ra£(f)), one 
can show that the solution of Eqs. and © 
can be written as, 


(n z (t)) =(n z ( 0)), (nl(t)) = 0, 

{m(t)) =e int $(t)(m( 0)), 
e int d$(t) 


<«^(*)> = - 


v cos 0 dt 


(m(0)>. 


( 6 ) 

( 7 ) 

( 8 ) 


Here we denote by $(t) the generating func¬ 
tional of the RTP 0 . For the FI decay, it is 

given by IEH3, 


<I> 7 (t) = e 74 sinh(7/rt) + cosh(7pf)^ , ( 9 ) 


where, p = \/l — {v cos d/7) 2 . 
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FI decay. We call noise weak if the dimen¬ 
sionless parameter, 77 = |tjcos 0 / 7 | -Cl. As it 
follows from Eq. ©> for weak noise, the de¬ 
cay rate of the non-diagonal averaged density 
matrix element (which characterizes the deco¬ 
herence) coincides with T v , provided by BR- 
theory. We call noise strong if 77 > 1. In 
particular, when 77 > 1 , the parameter, /i in 
Eq. © becomes imaginary, and the decay of 
the functional, <I>^(£), is accompanied by oscil¬ 
lations with frequency, 7 |/r|. 

Below, we compare the analytical solutions 
© ©, when the transverse effective field (re¬ 

laxation) is neglected, with the corresponding 
exact solutions obtained numerically from Eqs. 
© and ©. In numerical simulations, we put 
h = 1. All energy-dimensional parameters are 
measured in ps -1 (lps -1 ~ 0.66meV), and time 
is measured in ps. 

In Fig. m the strongly coupled dimer was 
considered (9 ss 0.968, e « 1.45). One can see, 
that in spite of the noise is strong: 77 ss 1.13 
(red dashed curve) and 77 Ri 2.27 (orange dashed 
curve), the approximate analytical solutions 
(shown by blue curves), are in a good agree¬ 
ment with the exact numerical solutions. In the 
inset, the noise amplitude is relatively strong 
(u = 20 ), and the matrix element, V 12 , of the 
Chll and Chl2 interaction, is large: e ss 8.24 
(strongly coupled dimer). At the same time, 
the noise is weak: 77 ss 0.24, and the BR ap¬ 
proach works. However, one cannot neglect the 
contribution from the transverse field to the de¬ 
coherence rate, = T v> + Yi/2. Indeed, in this 
case, T v « 0.3 and Ti « 0.48. That is why the 
approximate solution in the inset (blue curve) 
deviates significantly from the exact numerical 
solution (red dashed curve). 

Our numerical simulations show that the an¬ 
alytical solution, given by Eqs. © ®, is in 

a good agreement with the exact numerical so¬ 
lution, up to the value of e « 1.72 (9 < 7 r/ 3 ), in 
spite of the dimer is already strongly coupled. 

Echo signals. For simplicity, we assume that 
the 7t/2 and all spin echo 7 r-pulses act practi¬ 
cally instantaneously. The spin echo pulse, ap¬ 
plied at the time r, rotates the wave function 
around the x-axis, by the angle ir. The corre¬ 
sponding analytical solution can be written as, 


3> e (f) 


$^(t), 0 < t < r, 
t > T, 


( 10 ) 
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FIG. 1: (Color online) Time dependence (in ps) 
of 77 i = \(m(t))\, for the FI signal. Blue curves: 
analytical results. Dashed curves: exact solution. 
Parameters: fl = 127, 9 = 0.968, 7 = 10, v = 20 
(red curve), v = 40 (orange curve). Inset: 9 = 1.45, 
v = 20 . 


where, 

3>g(£) =&(t) + e _7t (l-(cosh( 7 /r(f - 2 r)) 

— cosh( 7 /if)). ( 11 ) 



FIG. 2: (Color online) Time dependence (in ps) of 
m = |(m(t))|. The FI decay: blue curve. Echo 
signals: green curve: N = 10; black curve: IV = 15 
(r = 0.025); cyan curve: N = 1; orange curve: N = 
10; red curve: N = 15 (r = 0.05). Parameters: 
f! = 10, v = 10, 7 = 10, 9 = 0. Inset: the sequence 
of echo 7 r-pulses applied at times: t = r, 2 r, 4r,.... 
The duration of the echo 7 r-pulse is: S = lfs, and 
its height is: h = 10 3 7r. 

In Fig. [2j the FI decay is shown for the trans¬ 
verse Bloch vector (blue curve). The N echo 
pulses are applied, partly restoring the FI sig¬ 
nal, and reducing the homogeneous broadening. 
Green curve: N = 10; black curve: N = 15 



























4 



0-1, , l . i 

0 0.5 1 1.5 


t 

FIG. 3: (Color online) Time dependence (in ps) 
of m = |(m(f))|. The FI decay: blue curve. Echo 
signals: red curve (number of echo pulses, N = 20). 
Parameters: 12 = 150, v = 40, 7 = 10, 6 = 0.165, 
t = 0.02 ps. Inset: 11 = 127, v = 20, 7 = 10, 
6 = 0.968, T = 0.03. 


(r = 0.025); cyan curve: N = 1; orange curve: 
N = 10; red curve: N = 15 (r = 0.05). Pa¬ 
rameters: 12 = 10, v = 10, 7 = 10, 0 = 0. The 
inset corresponds the the sequence of echo n- 
pulses applied at times: t = r, 2r, 4r,.... The 
duration of each echo pulse is: S = lfs, and its 
height is: h = 10 3 7 t . In this case, the noise is 
strong: r) = 1, and the dimer is weakly cou¬ 
pled: e = 0. (There is no contribution to the 
dynamics from the transverse field.) 

In Fig. [3] the time dependence (in ps) of 
the transverse Bloch vector, m = |(m(t))|, is 
shown. The FI decay corresponds to the blue 
curve. Echo signals: red curve (number of echo 
pulses, N = 20). Parameters: 12 = 150, v = 40, 
7 = 10, 9 = 0.165, r = 0.02. This corresponds 
to a weakly coupled dimer (e « 0.17), and to 
strong noise: rj m 3.9. The FI decay exhibits 
oscillations with the period: T = 27t/'T'ImI ~ 
0.16. The inset corresponds to the choice of 
parameters: fl = 127, v = 20, 7 = 10, 9 = 
0.968, r = 0.03. The noise is strong in this case, 
rj « 1.13, and the dimer is strongly coupled, 
e « 1.45. As one can see from Figs. [2] and [3l 
the homogeneous broadening of the FI signal 
can be significantly improved by applying the 
spin echo pulses for strong noise and for both 
weakly and strongly coupled dimers. 

Simultaneous action of the homogeneous and 
inhomogeneous disorder. Here, in addition to 
dynamical fluctuations (noise, £(2)), we con¬ 
sider an ensemble of TLSs (dimers) with fluc¬ 
tuating parameters, (ST, 6*, (f), due to the static 


disorder. It is well-known that this leads to 
the inhomogeneous broadening of the FI sig¬ 
nal decay. In our numerical simulations we as¬ 
sumed the independent Gaussian disorder for 
parameters, (12Our numerical simula¬ 
tions demonstrate that the main contribution 
from the static disorder, for a wide range of pa¬ 
rameters, is due to the fluctuations of the fre¬ 
quency, 12. So, below we neglect the static fluc¬ 
tuations of both angles, 9 and ip. We assume a 
Gaussian distribution for the random parame¬ 
ter 12, denoting the dispersion by er. Note that 
the results can easily be extended by including 
the static fluctuations of angles, 9 and ip, in the 
numerical solutions of the exact Eqs. a and 
©. 

In Fig. [4] we compare the results for the 
decay of the FI signal for three values of the 
dispersion, a = 0 , 10 , 20 , of static fluctuations 
of the frequency, 12 , and for the amplitudes of 
noise, v = 20,40. The parameter, 77 ss 1.13, 
so the noise is strong in this case. In the in¬ 
set, oscillations are observed with the period, 
T = 27r/7|/i|, due to the imaginary value of the 
parameter, p. 



FIG. 4: (Color online) Simultaneous action of ho¬ 
mogeneous and inhomogeneous broadening on the 
FI decay. The time dependence (in ps) of the gen¬ 
erating functional, 4>(f). Red curve: 0 = 0. Blue 
curve: a = 10. Green curve: a = 20. Parameters: 
12 = 127, 9 = 0.968, 7 = 10, v = 20. Inset: v = 40. 


The analytical solution, which includes the 
contributions from both, the static disorder and 
the dynamical noise, and corresponds to the 
spin echo signal applied at the time r, can be 
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written as, {m(t)) = $ < t(t)(m( 0)), where 


2+2 


a z t 


K(t) = 


2 <10 < t < t 

a 2 (t — 2 t ) 2 

t>T. 


( 12 ) 



FIG. 5: (Color online) The decay of the FI signal 
in the presence of both homogeneous and inhomo¬ 
geneous broadening (blue curve), and the action 
of the echo signal (red curve). Time dependence 
(in ps) of m = |(m(i))|. Parameters: = 127, 

6 = 0.968, 7 = 10, v = 20, a = 10, r = 0.075. 

In Fig. [U both static disorder of f2 and the 
dynamical noise, £(f), are included. The decay 
of the FI signal is shown by the blue curve. For 
our chosen parameters, the FI signal decays in 
approximately 200fs. The spin echo signal was 
applied at r = 75fs. As one can see, the echo 
pulse restores significantly the FI decay (red 
curve). Note that for the parameters chosen 
in Fig. [5j both dimensionless decay factors 
coincide at the characteristic time of the FI 
decay, f* = 200fs: yt* = cr 2 t*/2 = 2. So, both 
homogeneous and inhomogeneous broadening 
are partly compensated in this case by the spin 
echo signal. 

We presented the analytical and numerical 
results for spin echo pulses for the two-level 
systems (TLSs) chlorophyll-based dimers in 
bio-complexes, embedded in noisy protein en¬ 
vironment. We have shown that even strong 
dynamical broadening can be suppressed signif¬ 
icantly by the spin echo pulses. This is impor¬ 


tant for many bio-applications at ambient con¬ 
ditions. We also demonstrated the restoration 
of the free induction decay signal by the spin 
echo pulses when both homogeneous and in¬ 
homogeneous broadening equally contribute to 
the free induction decay. The application of the 
spin echo technique is especially useful for bio¬ 
systems with strong low-frequency noise, such 
as 1// noise. Our approach can be easily gen¬ 
eralized for this case, as was done in [7(, by 
introducing the corresponding ensemble of the 
fluctuators. Also, many different sources of dy¬ 
namical noises can be included in the presented 
here approach (as was done in 0 )- 
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